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Abstract 

We have extended our published set of low mass AGB stellar models to lower 
metallicity. Different mass loss rates have been explored. Interpolation formulae 
for luminosity, effective temperature, core mass, mass of dredge up material and 
maximum temperature in the convective zone generated by thermal pulses are 
provided. Finally, we discuss the modifications of these quantities as obtained 
when an appropriate treatment of the inward propagation of the convective in- 
stability, caused by the steep rise of the radiative opacity occurring when the 
convective envelope penetrates the H-depleted region, is taken into account. 
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1 Introduction. 

Results from detailed AGB stellar models demonstrate that low mass AGB stars (i.e. 
1.3 < M/Mq < 3) are the producers of the main component of the cosmic s-elements 
(Straniero et al. 1995; Gallino et al. 1998; Busso et al. 1999). Such a scenario has been 
confirmed by the measurements of the chemical composition of AGB stars (Lambert et 
al. 1995; Busso et al. 2001; Abia et al. 2001,2002) and by the analysis of the isotopic 
composition of meteoritic SiC grains (Gallino et al. 1997). 

The present generation of AGB stars, as observed in the disk of our Galaxy, have 
a nearly solar chemical composition. However, SiC grains found in pristine meteorites 
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originated in the C-rich circumstellar envelope of a pre-solar generation of AGB stars, 
whose original metallicity could have been somewhat lower than that found in the solar 
system material. In addition, an important contribution to our understanding of AGB 
stars comes from observations of the stellar population in the fields of the Small and 
the Large Magellanic Clouds, whose average metallicities are about 1/5 and 1/2 of the 
solar, respectively. In a previous paper (Straniero et al. 1997) the properties of low 
mass AGB stellar models with solar composition have been extensively discussed. Here 
we present an extension of these models to lower metallicity. 

2 The grid of models and interpolations. 

The full grid of models is summarized in Tab. 1. The corresponding AGB evolutionary 
sequences have been obtained by means of the FRANEC code, the same described in 
Straniero et al. (1997). The initial mass, the helium content and the metallicity are 
reported in column 1, 2 and 3, respectively. Various mass loss rates have been applied 
during the AGB phase. In all cases, we have used the Reimers formula (Reimers, 1975): 

M(M Q /yr) = 1.34 • 1(T 5 • V ■ j^- 

eff 

where M and L are in solar units and rj is a free parameter (see fourth column in 
Tab. 1). We have not considered mass loss in pre- AGB evolution, since, in the mass 
range 1.5-3 M , just a few hundredths of solar masses are expected to be lost. For 
example, by adopting r/=0.4 as representative of the pre-AGB mass loss rate (e.g. Fusi 
Pecci & Renzini 1976), we found that the total mass lost before the onset of the AGB 
phase by a 2 M star with solar composition should be 0.04 M Q . Clearly, such a small 
amount of mass loss has negligible consequences on the pre-AGB evolution. 

The thermally pulsing AGB models (TP-AGB) are characterized by three distinct 
phases: 

1) The early phase. It includes the first few (gentle) thermal pulses which are not 
followed by the third dredge up (TDU); 

2) The TDU phase. It begins when the mass of the H exhausted core becomes as 
large as 0.59 — 0.6M Q . The amount of mass dredged up firstly increases, as the core 
mass increases, and then decreases, when the envelope mass is substantially reduced; 

3) The final phase. During the last part of the AGB, the concurrent actions of 
mass loss and H-burning erode the envelope mass and the TDU ceases. The minimum 
envelope mass for the occurrence of the TDU is about 0.4 — O.5M . 

The variation of the mass dredged up, as a function of the core mass, is shown in 
Fig. 1 for some selected sequences of models. 

An interpolation on the grid of models provides useful analytic expressions for 
relevant quantities, as a function of M H (the mass of the H-exhausted core, in M ), 
M env (the mass of the H-rich envelope, in M ) and Z (the initial mass fraction of 
elements with A > 12). 

In the following we report the result of these interpolations. 
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Figure 1: The amount of mass dredged up per thermal pulse for three evolutionary 
sequences: M = 2.0 M , Z = 0.02 and r) = 0.4 - triangles (FRANEC models) and solid 
line (interpolation formula); M = 1.5 M , Z = 0.006 and rj = 0.5 - circles (FRANEC 
models) and dashed line (interpolation formula); M = 1.5 M , Z = 0.003 and rj = 0.4 
- squares (FRANEC models) and dotted line (interpolation formula). 
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1. The time elapsed between two successive thermal pulses (interpulse period): 
At ip (10 5 yr) = 18.5178 - 59.8876 • M H + 48.8462 • M 2 - 4.0273 • \ogZ+ 

+5.8422 • M H • log Z 

which is valid for Mh > 0.58. For smaller core masses, substitute Mh with (1.16 — 
Mb). The standard deviation of this fit (i.e. the average difference between the values 
calculated with the fitting formula and those of the models) is 0.01. 

2. The mass of the H-depleted material dredged up in a given TDU episode: 
AM TDU (M G ) = (1.0 + 0.21 • log ^ + 6.3 • (log ^) 2 )(-6.2 • 10- 4 + 

+9.1 • 10- 4 • M env - 3.7 • 10" 4 • M 2 env + 5.061 • 10~ 2 • M env ■ 5M n - 

-5.96 • 10~ 3 • M 2 env ■ 5M U - 3.8372 • 10' 1 • M env ■ 5M n + 9.448 • 10" 2 • M 2 env • <5M 2 ) 

where Z Q = 0.02 and <5Mh = — M u and M H is the core mass at the last thermal 
pulse without TDU (last pulse of the early phase) . For the masses and chemical com- 
positions we have considered, ranges between 0.59 and 0.6. The curves reported 
in Fig. 1 have been obtained by using = 0.595. When M env < 0.4 or 5M U < 0, 
take AM TDU = 0. The standard deviation of this fitting formula is 2 • 10" 4 . 

3) The luminosity at half of the interpulse period: 

log ^ = 3.603 + 6.7806 • AM H - 27.582 • AM 2 + 277.333 • AM 4 

here AM H = M H — M H (to) and M H (to) is the core mass at the epoch of the first 
thermal pulse (to). The latter quantity slightly depends on the metallicity: M H (to) = 
0.517 - 1.93910" 2 x logZ. For this fit the standard deviation is 10~ 4 . 

4) The effective temperature at half of the interpulse period: 

\ogT eff (K) = 5.0475+0.1438dogZ e// +0.0513-(logZ e// ) 2 -4.1895-M H +2.9594-M 2 

here, Z e ff is the effective metallicity, which includes the extra-carbon in the envelope 
due to the third dredge up. The use of Z e ff, instead of Z, has minor effects on solar 
composition stars, but may significantly increase the estimated radius and, in turn, the 
mass loss rate for stars with an original low content of metal. The standard deviation 
is 5 • 10~ 4 . 

5) The maximum temperature attained at the bottom of the convective zone generated 
by a thermal pulse: 

\ogT max (K) = 6.7747 + 4.6856 • M H - 3.21 • M 2 + 0.01 • M env + 

-6.2337 • Z + 7.87595 • Z ■ M H 
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for \ogT max the standard deviation is 3 • 1CT 3 . Finally, the evolution of the core 
mass during the whole TP-AGB phase may be estimated by means of the following 
formula: 

Mh = My + 6.9 • KT 3 • k - 6.0 ■ 10~ 5 • k 2 

where k = —n,—n+ 1, 0, 1, j. k = corresponds to the last thermal pulse 
without TDU (Mh = M H ); n+ 1 is the number of TPs of the early phase, while j is the 
number of TPs occurred from the first TDU episode up to the AGB tip (TDU phase 
+ final phase). This formula gives M H within 5 -1CT 3 M . Then, the envelope mass is 
easily obtained: M env = M — M H , where M is the total mass. 

Note that, if the initial mass of the star is too small, the conditions for the activation 
of the third dredge up (namely, Mh > 0.60 M and M env > 0.4 M Q ) are never reached. 
The minimum mass for the occurrence of the TDU, as a function of the mass loss rate, 
are reported in Fig. 2, for solar composition models. In stars with masses exceeding 
this lower limit, the envelope composition is modified by the TDU: 12 C and s-elements 
are enhanced. Obviously, a C-star may form, after a suitable number of TDU episodes, 
only if the initial mass is large enough. The dashed line in Fig. 2 shows, for solar 
metallicity models, the lower mass limit for the formation of a C-star as a function of 
V- 

Both these lower limits depend on metallicity: the smaller the metallicity, the 
smaller is the mass for the activation of the TDU phase and the lower is the minimum 
initial mass for the formation of a C-star. This occurrence is illustrated in Fig. 3 (for 
77=0.5). 

At low metallicity (i.e. Z < 4 • 10~ 3 ), even taking into account an enhanchement of 
the a-elements (as typically found in halo stars), the amount of oxygen in the envelope 
is so low that just one TDU episode is sufficient to increase the C/O ratio above unity. 
Correspondingly, the core mass and the luminosity at C/O = 1 decrease at low Z (see 
Fig. 4). 

3 The efficiency of the third dredge up: the Opacity 
Induced Convection. 

An He-rich zone, as the one left by the shell H-burning, has a significantly lower 
opacity than an H-rich zone, at the same temperature and pressure. Thus, when 
the innermost layer of the convective envelope (H-rich) penetrates into the region of 
decreasing H and increasing He (at the epoch of the TDU), the local value of the 
radiative opacity increases. Since the radiative gradient (V mc z) is proportional to the 
opacity, the condition V ra d ^> V ' a d (where V ' a d is the adiabatic temperature gradient) 
takes place at the base of the convective envelope. In this case, the internal boundary 
of the convective zone becomes unstable (Becker & Iben 1979; Castellani, Chieffi & 
Straniero 1990; Frost & Lattanzio 1996; Castellani, Marconi & Straniero 1998). In fact, 
even a small perturbation may increase the radiative gradient in the formally stable 
region, located immediately below the convective envelope, which immediately becomes 
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Figure 2: The minimum initial mass for the occurrence of the TDU (solid lines) and 
the minimum initial mass for C-star progenitors (dashed lines) as a function of the 
mass loss rate. The 77 values reported in the abscissa refer to the free parameter in the 
Reimers formula. These results are based on solar metallicity models (Z = 0.02 and 
Y = 0.28). Thick lines represent the case of models with minimal dredge up, while 
thin lines have been derived from models that include an exponential decline of the 
convective velocity at the boundaries of the convective zone (see section 3). 
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Figure 3: The minimum initial mass for the occurrence of the TDU (solid line) and the 
minimum initial mass for C-star progenitors (dashed line) as a function of Z. These 
results are based on models computed with a Reimers mass loss rate (rj = 0.5) and 
minimal dredge up (see section 3). 
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Table 1: The computed grid of models. M is the initial mass. Y and Z are the 
initial helium content and the metallicity, respectively. r\ is the free parameter in the 
Reimers mass loss formula. (5 is the free parameter in the exponential decline law for 
the boundary velocity of the convective elements (section 3). 



convectively unstable. At present, a satisfactory treatment of this phenomenon has 
not yet been found. However, its overall effect on AGB stellar models may be derived 
by appling a small perturbation to the boundaries of the convective regions. As a 
tentative scheme, we have supposed that the velocity of the convective elements, usually 
evaluated by means of the mixing length theory, does not abruptly drop to at the 
convective boundaries, but decreases following an exponential decline law (for more 
details on the mixing algorithm see section 2 in Chieffi et al. 2001): 

v = v bce exp (-^), 

where z is the distance from the convective boundary, t> bce the average velocity at the 
convective boundary, H P the pressure scale height at the convective boundary and (3 is 
a free parameter that controls the steepness of the velocity decline. Usually, when the 
convective boundary is located within a chemically homogeneous region, vj, ce ~ and 
the convective border is stable. However, at the epoch of the TDU, at the base of the 
convective envelope the difference between radiative and adiabatic gradients increases. 
Thus, Vbce increases, causing an inward propagation of the convective instability. The 
final result is a substantial increase of the TDU efficiency. With (3 = 0.1, for example, 
the amount of H-depleted mass dredged up by convection is about a factor 2 larger 
than that found in models without the exponential decline of the convective velocity 
(Cristallo et al. 2003). 

Few additional models have been computed by taken j3 — 0.1 (see last column in 
Tab. 1). The effect on the minimum mass for C-star progenitors is shown in Fig. 2, 
in the case of solar composition models. The corresponding decrease of the minimum 
C-star luminosity is illustrated in Fig. 4. 
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Figure 4: The minimum luminosity for a C-star as a function of Z. These results are 
based on models computed with a Reimers mass loss rate (rj = 0.5) and with (dashed 
line) and without (solid line) exponential decline of the convective velocity. 



4 Final remarks 



In this paper, based on an homogeneous grid of low mass AGB stellar models, we 
provide some interpolation formulae for basic ingredients used in nucleosynthesis cal- 
culations. These formulae are strictly valid for M < 3M , at solar metallicity, and for 
M < 2.5M , in the metallicity range 0.003 < Z < 0.02. 

A study of the third dredge-up as a function of total mass and metallicity has 
been recently published by Karakas, Lattanzio & Pols (2002), while a comparison 
among the AGB stellar models computed by means of different evolutionary codes has 
been discussed by Lugaro et al. (2003). Once the same mass loss rate is adopted, 
a substantial agreement is found between our (FRANEC) results and those obtained 
by using the Mount Stromlo Stellar Structure Program (Frost and Lattanzio, 1996; 
Karakas, Lattanzio & Pols 2002). 

As a final remark, note that although we have used a Reimers formula to evaluate 
the mass loss, since the fittings formulae are given as a function of the envelope and 
core masses, they can be used for any mass loss law. 
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